function y = initial_function2(x)
y = sin(x) + sin(100*x);